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J^T ■ Abstract 
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The statistical coalescence model for the production of open and hidden charm 
is considered within the canonical ensemble formulation. The data for the J ftp 
multiplicity in Pb+Pb collisions at 158 A-GeV are used for the model prediction of 
the open charm yield which has not yet been measured in these reactions. 

The charmonium states J/ip and ip' have been measured in nucleus- nucleus (A+A) 
collisions at CERN SPS over the last 15 years by the NA38 and NA50 Collaborations. 
This experimental program was motivated by a suggestion [jXJ] to use the J/tp as a probe of 
the state of matter created in the early stage of the collision. In this approach a significant 
suppression of J / ip production relative to Drell-Yan lepton pairs is predicted when going 
from peripheral to central Pb+Pb interactions at 158 A-GeV. This is originally attributed 
to the formation of a quark-gluon plasma, but could be also explained in microscopic 
hadron models as secondary collision effects (see 0] and references therein). 

The statistical approach, formulated in Ref.||, assumes that J/ip mesons are created 
at hadronization according to the available hadronic phase-space. In this model the J/ip 
yield is independent of the open charm yield. The model offers a natural explanation of 
the proportionality of the J/ip and pion yields and the magnitude of the J/ijj multiplicity 
in hadronic and nuclear collisions. 

Recently the statistical coalescence model M and the microscopical coalescence model 
|| were introduced for the charmonium production. Similar to the statistical model ||, 
the charmonium states are assumed to be formed at the hadronization stage. However, 
they are produced as a coalescence of created earlier c-c quarks and therefore the multi- 
plicities of open and hidden charm hadrons are connected in the coalescence models. In 
Ref. dJ the charm quark-antiquark pairs are assumed to be created at the early stage of 
A+A collision and the average number of c-c pairs, N^ r , is fixed by the model consid- 
eration based on the hard scattering approach. The estimated number iV* r seems to be 
larger than the equilibrium hadron gas (HG) result. This requires the introduction of a 
new parameter in the HG approach [|J - the charm enhancement factor 7 C (it was denoted 
as g c in Ref.H). This is analogous to the introduction of strangeness suppression factor 
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7s in the HG model, where the total strangeness observed is smaller than its thermal 
equilibrium value. Within this approach the open charm hadron yield is enhanced by a 
factor 7 C and charmonium yield by a factor 7^ in comparison with the equilibrium HG 
predictions. 

The thermal HG calculations in Ref . Q are done within the grand canonical ensemble 
(g.c.e.) formulation although the validity of the g.c.e. results for the charm hadron yield 
was questioned by the authors of [f| themselves. As the total number of charm hadrons 
is expected to be smaller than unity even for the most central Pb+Pb collisions an exact 
charm conservation within the canonical ensemble (c.e.) should be imposed^. Note also 
that the c.e. formulation was successfully used in Ref. JF) to calculate the open charm 
hadron abundances in e + e~ collisions at ^/s = 91.2 GeV with an experimental input 
of the total open charm production. In this letter stimulated by the above proposals 
we consider the c.e. HG formulation for the physical picture suggested in Ref.[[|]. The 
experimental value of the J/ip multiplicity (J/ip) will be then used to predict the open 
charm yield within the statistical coalescence model. 

The main assumption of Ref. 0] is formulated as 

NM r = \lcN Q + l 2 c N H , (1) 

where No is the total thermal multiplicity of all open charm and anticharm mesons and 
(anti)baryons and Nh is the total thermal multiplicity of particles with hidden charm. 
Note that open charm resonance states (not included in |4[]) give essential contribution^ to 
Nq- The number of directly produced cc pairs N^ r in the hard collisions is estimated in 
Ref.§ to be equal to iVj r = 0.17 for Pb+Pb SPS collisions with N p = 400 participants. 
This number is however not quite confident. The average number of cc pairs created 
in nucleon-nucleon collisions at = 17.3 GeV was estimated from existing data as 
3 • 10~ 4 in Ref. [§]. With a linear on N p extrapolation to central A+A collisions one 
obtains N% r = 0.06 for N p = 400, but assuming that open charm production in central 
A+A collisions scales as N^ 3 an estimate N^ r 0.35 is obtained for N p = 400 [|. Note 
also that a recent analysis of the dimuon spectrum measured in central Pb+Pb collisions 
at 158 A-GeV by NA50 Collaboration || suggests a significant enhancement of dilepton 
production in the intermediate mass region (1.5-^2.5 GeV) over the standard sources. 
The primary interpretation attributes this observation to the enhanced production of 
open charm ||: about 3 times above the pQCD prediction for the open charm yield in 
Pb+Pb collisions at SPS. 

In Ref. [U No and Nh are calculated in the g.c.e.. In the g.c.e. the thermal multiplic- 
ities of both open charm and charmonium states are given as (Bose and Fermi effects are 
negligible): 

djV e^ /T m „ T ^ /toA , Tr „ /T /m,T\ 3/2 / m,\ 

Nj = 2** T m? K2 vt) dj v 6 vir) exp v~t) ' (2) 

where V and T correspond to the volume^ and temperature of HG system, rrij, dj denote 

1 This was first suggested by K. Redlich and B. Miiller (quoted in Ref.Q and L. McLerran, private 
communication) . 

2 We are thankful to Braun-Munzinger and Stachel for pointing this out. 

3 To avoid further complications we use ideal HG formulae and neglect excluded volume corrections. 
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particle masses and degeneracy factors and Ki is the modified Bessel function. The 
particle chemical potential fij in Eq.(fJ) is defined as 

N = b jf i B + Sjfi S + Cjfi C , (3) 

where bj , Sj , Cj denote the baryonic number strangeness and charm of particle j . The 
baryonic chemical potential /i# regulates the baryonic density of the HG system whereas 
strange us and charm fic chemical potentials should be found from the requirement of 
zero value for the total strangeness and charm in the system (in our consideration we 
neglect small effects of a non-zero electrical chemical potential). 

In the c.e. formulation (i.e. when the requirement of zero "charm charge" of the HG is 
used in the exact form) the thermal charmonium multiplicities are still given by Eq. (Q) as 
charmonium states have zero charm charge. The multiplicities (Q) of open charm hadrons 
will however be multiplied by an additional 'canonical suppression' factor (see e.g. [|T0] ). 
This suppression factor is the same for all individual open charm states. It leads to the 
total open charm multiplicity Nq in the c.e.: 

N ce ^_ at hjNp) 

where No is the total g.c.e. multiplicity of all open charm and anticharm mesons and 
(anti)baryons calculated with Eq.(|2]) and Iq, I\ are the modified Bessel functions. For large 
open charm multiplicity No » 1 one finds Ii(No) / 'Iq(No) — * 1 and therefore Nq — > No, 
i.e. the g.c.e. and c.e. results coincide. For Nq « 1 one has Ii(Nq) / Iq(Nq) — Nq/2 
and Nq = Nq ■ Nq/2, therefore, Nq is strongly suppressed in comparison to the g.c.e. 
result N ■ 

Assuming the presence of the charm enhancement factor 7 C the statistical coalescence 
model within the c.e. should now be formulated as: 

JV? = \ 7o No + N H . (5) 

Therefore, the baryonic number, strangeness and electric charge of the HG system are 
treated in our approach according to the g.c.e. but charm is considered in the c.e. for- 
mulation where the exact charge conservation is imposed. 

The logic of Ref.|4j] is the following: 1) input N^ r number (it is assumed to be equal 
to 0.17 for N p = 400) into Eq.(|l]); 2) calculate the 7 C value; 3) obtain J/ip multiplicity 
as (J/?p) = j^Nj/^, where Nj/^p is given by Eq.(|). Note that the second term in both 
Eq.(|IJ) and (||) gives only a tiny correction to the first term. Therefore, 7 C = 2N^ r /No- 

Our consideration differs from that of Ref. Q in three points. First, we use Eq.@ 
instead of ([]]). Second, in our calculations we take into account all known particles and 



resonances with open and hidden charm [jTTJ. Third, we will proceed with Eq.(|5|) in the 



reverse way. As the (J ftp) multiplicities can be extracted from the NA50 data on Pb+Pb 
collisions at 158 A-GeV for different values of N p , we start from the requirement: 

(m = ii n i j% , (6) 
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to fix the 7 C factor. In Eq.(^|) the total J/ip thermal multiplicity is calculated as 

Nj% = Nj/^ + BrW)Nv + Br( X i)N Xl + Br( X2 )N X2 , (7) 

where Nj/^, N f >, N Xl , N X2 are given by Eq.(@) and Br(i//) = 0.54, Br{ X i) = 0.27, 
Br(x2) — 0.14 are the decay branching ratios of the excited charmonium states into J/if>. 
Eq.(H) will be used then to calculate the value of iVjT. This value will be considered as 
a prediction of the statistical coalescence model: the open charm yield has not yet been 
measured in Pb+Pb collisions at SPS. 

We use two different sets of the chemical freeze-out parameters: 



A: T = 168 MeV, y, B = 266 MeV, 7, = 1 Ref.pf; 

B: T = 175 MeV, fi B = 240 MeV, 7, = 0.9 Ref. || . 

They both were fixed by the HG model fit to the hadron yields data in Pb+Pb collisions 
at 158 A-GeV (the inclusion of open charm and charmonium states does not modify the 
rest of the hadron yields). For the fixed number of participants N p the volume V is defined 
then from N p = Vtib, where n B = n B (T, ^Bils) is the baryonic density calculated in the 
g.c.e.. With two sets of the chemical freeze-out parameters A and B we find Nj°L and Nq 
values using Eq.(|J), calculate 7 C factors from Eq.(|6|) and then calculate Nf~ r from Eq.(H). 

In the above c.e. consideration with exact charm conservation the 7 C parameter reg- 
ulates the average number of cc-pairs in the HG. Therefore, N c = is restricted 
exactly (the c.e.), but the value of N c + iV E is restricted on average (the g.c.e). The above 
c.e. calculations are based therefore on the thermal model distribution for probabilities 
to observe 0, 1, 2, ... of cc-pairs in the equilibrium HG. One needs then an additional pa- 
rameter 7 C to adjust these thermal probabilities to the required number of N^ r . Another 
way is to restrict also the N c + N% numbers in the c.e. calculations and use non-thermal 
probabilities to create k — 1,2, ... of cc-pairs in hard collisions. For the fixed number of 
cc-pairs equal to k, the average multiplicity of hidden charm can be approximately calcu- 
lated in the following way. We keep in the c.e. HG partition function the leading terms 
only with and 1 hidden charm particles and neglect all configurations with 2,3,..,fc char- 
monium particles. This corresponds to the expansion in powers of the small parameter 
N H /(N /2) 2 « 1. It gives: 

(N H ) k ~ e Nh 



(No/2) 2 ' 

where multiplicities N Q and N H are calculated in the g.c.e. using Eq.(Q). 

Because of the assumed hard scattering origin of the cc production the Poisson dis- 
tribution P(k) = f k exp(—f)/k\ looks quite natural (k = 0, 1, 2, ... is the number of pairs 
created, / = iVjT is the average number of pairs). The calculations with these 'dynamical' 
probabilities contain no additional free parameter. All 'dynamical' information needed 
for the c.e. calculation is now given by the value of / (parameter 7 C does not appear). 
With Eq.(||) the result for J/ijj yield is: 

]\f tot 
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where Nj°K is given by Eq.(|7]). Using the experimental values for (J/if>) one obtains from 
Eq.(^|) the average number of cc-pairs / = iVjT. The results appear to be rather close 
to those obtained with thermal probabilities. The reason of this fact is that states with 
k = and k — 1 dominate in both thermal and 'dynamical' probability distributions. To 
illustrate this lets consider an extreme choice: the HG states with N c = = 1 appear 
with probability / = JVjj r , the HG states with N c = = appear with probability 
1 — /, and states with more than one cc-pairs are neglected. Under these restrictions the 
J /if) multiplicity becomes equal to: 

< w > = ! wfe ■ (10) 



One sees that Eq.flTUD is close to Eq.(|) if / « 1 and N H « (N /2) 2 . From Eq.(|lOD 
one finds: 



/ = J^T (No/2) 2 + N H = Y c (N /2Y + N H 



(JM 

Ttot 

This coincides with Eq.(j5l) at small values of JcNq- 



'ID 



We present now the model calculations for central Pb+Pb interactions at 158 A-GeV 
(N p = 100 -7- 400). Using the estimates for experimental J /if) multiplicities and assuming 
that the system volume V scales linearly with N p (i.e. N p = Vub(T, [Ib, Is)) we ; first, 
calculate the thermal J /if) multiplicity N^h (|^) - including the feeding from the excited 
charmonium states. Then we use Eq.(|5]) to find the parameter j c . Finally, we calculate 
the predicted values of Nffi r from Eq.(|^). The results are presented in Tables 1 and 2 
where the sets of the chemical freeze-out parameters A and B are respectively used. 

A reliable extraction of the J /if) yields from the published data appears to be non- 
trivialQ, Ref.[14j] suggests an approximately linear increase of (J /if)) with N p . The results 



for (J /if)) presented in Ref. || were evaluated from the data of the NA50 Collaboration 
fTB| using the procedure described in These results for {J /if)) are used as input for 



the statistical coalescence model analysis in Tables 1 and 2. Assuming that iVjT scales 
as we find from Tables 1 and 2 a value of a — 1.6 -v- 1.7. This value is larger than 
a = 4/3 expected in the hard-collision model. Although the values of Njfy, No and 7 C 
are rather sensitive to the temperature parameter, the model predictions for Nf}/ r 
essentially unchanged. 



In conclusion, the statistical coalescence model with an exact charm conservation is 
formulated. The canonical ensemble suppression effects are important for the thermal 
open charm yield even at N p = 400. These effects become crucial when the number 
of participants N p decreases. From the J /if) multiplicity data in Pb+Pb collisions at 
158 A-GeV the open charm yield is predicted: N^ r = 0.5 -r- 0.6 in central (N p = 360) 
collisions. An uncertainty of this prediction is mainly because of uncertainties in different 
compilations of the (J/ if) data. 

4 We are thankful to K. Redlich and M. Gazdzicki for the useful comments. 
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Table 1 
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N l j% ■ 10 4 
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NA50 data 


Eq.® 
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Poisson 




Compil. § 


Set A 


Set A 


Eq.© 


Eq.® 


Eq.(|) 


100 


2.2±02 


0.56 


0.26 


2.0 


0.066 


0.064 


200 


3.9±0.2 


1.1 


0.52 


1.9 


0.21 


0.20 


300 


6.4±0.6 


1.7 


0.79 


2.0 


0.46 


0.41 


360 


6.9±0.7 


2.0 


0.94 


1.9 


0.57 


0.51 



Table 2 
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Poisson 
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Eq.® 


Eq.® 
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100 


2.2±02 


1.1 


0.39 


1.4 


0.072 


0.070 


200 


3.9±0.2 


2.2 


0.77 


1.3 


0.23 


0.22 


300 


6.4±0.6 


3.3 


1.17 


1.4 


0.50 


0.45 


360 


6.9±0.7 


4.0 


1.40 


1.3 


0.62 


0.55 



It is interesting to compare our estimate N^ r = 0.5 -~ 0.6 with results predicted in 
different model approaches. The pQCD inspired models suggest the values of N^ r = 
0.1 0.3 (the value of N^ r = 0.17 is an estimate of Ref.[|). Much larger value of 
N c5 = 3.4 is obtained in Ref. within the microscopic coalescence model. Even larger 



value of N cS ~ 8 is suggested in Ref. |L6 assuming the charm equilibration in the quark- 



gluon plasma at the very early stage of Pb+Pb reaction. 

The statistical coalescence model predicts also the N p dependence of N^ r and the 
yields of individual open charm states. All these predictions of the statistical coalescence 
model (the open charm yield has not been measured in Pb+Pb) can be tested in the 
near future (measurements of the open charm are planned at CERN). This will require 
to specify more accurately the (J/if>) data. 

The charm enhancement factor 7 C found from the (J/ip) data appears to be not much 
different from unity and its value is rather sensitive to the temperature parameter. There- 
fore, both the statistical model of Ref.|| and the statistical coalescence model considered 
in the present paper lead to similar results for the J /if} yield. However, the predictions of 
these two models will differ greatly at RHIC energies: according to @ the J/if> to pion 
ratio is expected to be approximately equal to its value at the SPS, but according to the 
statistical coalescence model this ratio should increase very strongly. The predictions of 
the present model for the RHIC energies will be presented elsewhere. 
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